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A transport theory which is not restricted to the gradient and quasi-particle 
Cn I approximations is presented which is formulated in terms of the energy moments, 

or equivalently the equal-time derivatives of the one-particle Green functions. A set 
(NJ , of infinite hierarchies of kinetic and constraint equations for equal-time quantities 

^^ I for the spectral and the kinetic part of the one-particle Green function is derived. 

f— ^ ' The hierarchies for the spectral function truncate automatically as in the mean field 

^P . approximation. The possibility of a systematic truncation of the hierarchies for 

the kinetic part is discussed. The effects of the quantum corrections are illustrated 

in a simple one-dimensional model. 
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1 Introduction 

A transport description of many-body non-equilibrium systems is a very convenient 
approximation used with success in the study of heavy ion colhsions in the last two 
decades [|1|. It has also been used for more than 40 years in solid state physics to 
describe the transport properties like electric conductivity, thermal conductivity, etc. In 
this approach the one-body phase-space distribution function satisfies a semi-classical 
transport equation of Boltzmann type. The validity of a transport description can be 
established by a rigorous derivation from quantum field theory 0. 

A main step in the derivation of transport equations starting from the Dyson 
equation for the Green function is to put particles on the mass-shell [§]. On the other 
hand, off-shell effects are thought to be important in heavy ion collisions for processes 
like subthreshold production of particles or for the inclusion of baryonic resonances in 
the transport description [Q. This is the case, in particular, for the A-resonance for 
which the spectral function has a large width even outside the medium. 

The off-shell effects are related to the dependence of one-particle Green functions on 
two time arguments. A sytematic way to include the off-shell effects into the dynamics 
of a particle in an external c-number field (or in the mean field approximation) was 
developed 0] resulting in transport and constraint equations for equal time quantities 
which are derivatives of the Green function in the relative time at zero relative time, 
or equivalently energy moments of its Wigner transform. The main advantage of this 
equal-time description is the possibility to specify the initial conditions necessary to solve 
the transport equations. 

More recently new derivations of semiclassical transport equations 0, ^, ^ do not 
use the on-shell approximation. They result in first order partial differential equation for 
some function defined on 8-dimensional (phase) space which is numerically solved using 
the test particle method by further approximating it by a sum of 6 functions. It has been 
argued in ref. that in general the function that can be represented by a test particle 
ansatz is not the phase space distribution function, but some more general quantity. 
This complicates the solution of the transport equations in 8-dimensional phase space, 
giving interest to alternative methods. In addition these treatments use the first-order 
gradient approximation which is based on the assumption of a weakly inhomogeneous 



and slowly evolving system. The gradient approximation may loose validity when the 
range of the interactions becomes comparable with that of the space inhomogeneities of 
the system. Other authors |§] therefore generalize the Boltzmann equation by including 
third order terms in the space gradients and obtain equations for quantities which depend 
on position, time, energy and two momentum components. 

In the present work we propose a systematic way to derive equations for quantities 
which characterize the off-shell behaviour of the Green functions and also include higher 
order gradient terms, i.e. go beyond the mean field approximation. The idea is to 
formulate a theory for the energy moments of the one-particle Green functions. These 
energy moments measure the shape of the Green functions in energy, e.g. the second 
order energy moment is a measure of the width of the spectral function. Generally the 
ra^'^-order energy moment is related to the n*^-order derivative of the Green function 
with respect to the relative time taken at zero relative time. We use the method of 
real-time Green functions defined on the Schwinger-Keldysh time path, which is the 
appropriate formalism to describe quantum systems in non-equilibrium. Taking the 
Wigner transform of the Dyson equation in space coordinates only we obtain evolution 
and constraint equations for the energy moments of Green functions. The result is a set 
of four infinite hierarchies of equations. A systematic way to truncate these hierarchies 
is discussed. 

As we are going beyond the first order gradient approximation, we can study the 
non-localization effects due to the quantum motion. There have been attempts to include 
the nonlocal terms in the scattering integral ^ and this leads to consistency with the 
thermodynamic virial corrections . 

In the last section we will illustrate in a simple 1-dimensional example the effects 
of the higher order gradient terms in the evolution equation as well as the influence of 
the memory effects on the velocity of thermalization. 



2 Spatial Wigner transforms of Dyson equations 

To keep the formalism simple we consider in this work a non-relativistic system of 
fermions described by the Hamiltonian 
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where '^\j and '^ h are the field operators in the Heisenberg picture and V is the two- 
particle interaction potential {% = !). A Dirac or Klein-Gordon field can be treated 
along the same lines, however, here we want to avoid the technical complications related 
to the Dirac matrix structure of the Green function and self-energy. 

For a non-equilibrium situation it is convenient to introduce a path-ordered Green 
function defined on the Schwinger-Keldysh double time contour |jll[]. The one-particle 
Green function on the contour G is defined as the expectation value of a product of field 
operators in the Heisenberg picture 
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where 1 = (xi,ti) and pu is the time-independent density operator which specifies 
the state. Depending on the position of the two time arguments on the contour G is 
equivalent to four different Green functions defined on the real time axis: causal and 
anticausal when the two time arguments are on the same branch of the contour and 
correlation Green functions when the time arguments are on different branches 
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The causal and anticausal Green functions satisfy the following decomposition relations 
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It is useful to define retarded and advanced Green functions 
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G+'-(l,l') = ±e{±{U-tr))[G>{l,l')-G<{l,l')] 
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as well as spectral and Keldysh (or kinetic) Green functions 

G^(l,l') = G>(1,1')-G'<(1,1') (6) 

G^(l,l') = G>(1,1')+G<(1,1'). (7) 

We note [Q] that all the Green functions defined above can be expressed in terms of only 
two independent Green functions, for which we take G^ and G^ . For future reference 
we note the relations 

(G^'^(l,l'))' = -G^'^(l',l) (8) 

and, as a result of the equal-time commutation relations, the property 

G^(l,l')|i,=i,, = 5(xi-xi0. (9) 

(Qc Q< \ 
„> „jj , and correspond- 
ingly for the other one-particle quantities, the Dyson equation can be written in the 
equivalent forms 

G(l,l') = Goil,l') + J dl' J dl"Goil,l")Ul",l"')Gil"',l') (10) 

G(l,l') = Go{l,l') + I dl' I dl"G{l,l")Ul",l"')Goil"',l') (11) 

where the time arguments take values on the Schwinger- Keldysh contour and Gq is the 
free Green function which satisfies the equations 

with 6_ the Dirac 5-function defined on the contour. The self-energy S contains all the 
information on the interaction. The way in which it is approximated specifies the type 
of approximation used. We will discuss this for a few cases below. Its causal and anti- 
causal components have the same decomposition relations as the Green functions, eq. 
(5), provided the singular, i.e. time-local, part on the contour is separated off [|T^ 

S(l,l') = S5(xi,xi,,ti)5(ti-tiO + S(l,l') (13) 



This singular part is given by the Hartree and Fock diagrams which are local in time, and 
are also discussed further in Sect. 5 below. Also for S we apply analogous definitions 
as in eqs. (6,7). In the following, to avoid writing integrals, we consider the space-time 
coordinates as matrix indices and use the Einstein convention to sum (integrate) over 
the repeated indices (coordinates). 

The equations satisfied by the interacting Green functions, eqs. (10) and (11), can 
be put in a triangular matrix form for G = ^^ „_,_ [|T^. From the non-diagonal 
part we obtain the equations for the kinetic part of the Green function 

siG'^l, 1') - 2^(1, l")G'\l", I') = S+(l, l")G'\l", I') + S^(l, ^^^-(l", 1') 

4G^(1, 1') - G'^'(l, 1")S5(1", 1') = G^(l, 1")S-(1", 1') + G+(l, 1")S^(1", 1') (14) 

where Sj = idt^ + -^^i- /,From the diagonal part we obtain analogous equations for the 
spectral part of the Green function, G^. 

To derive an equal-time formulation we perform a Wigner transform of these equa- 
tions in the space coordinates only. For the Wigner transform of a general function F 
we use the notation 
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The transforms of G and S are denoted by g and a, respectively. 
Introducing the differential operator (Poisson operator) 

<>{f-9) = (Vi?/)(Vp^)-(Vp/)(V/?^) 
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the equations for the Wigner transformed spectral and kinetic Green functions, g^ and 



g^ , are 
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where we did not indicate the dependence of g"^'^ and a"^'^ on R and p. The upper and 
lower signs in eqs. (17-19) resuh when starting from eqs. (10) and (11), respectively. 
The boundaries of the integrals in eqs. (18,19) come from the step function in eq. (5). 

Setting a^ = a^ = 0, i.e. when the scattering is neglected, the problem is treated 
in the mean field approximation, the terms on the right-hand side of the above equations 
vanish, and eqs. for g"^ and g^ are identical. Therefore one has only two equations for 
g^ = ^{g^ — g^) which is directly related to the particle density. In the general case 
it is neccessary to consider all the four equations. Taking the sum and the difference of 
eqs. (17) we obtain for g^ and g^ two sets of equations 
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The equations in the upper row are called kinetic equations because they involve the 
time evolution in the global time T. The ones in the lower row are called constraint 
equations [Q, since under certain conditions (see below) they are algebraic equations. 



3 Hierarchies of equations for the spectral functions 

The eqs. (20) are, of course, still exact. The quantities g^'^{T;T), a^'^(T;T) depend 
on two time arguments, the global time T and the relative time r. The integrations in 
the collision terms, (eqs. (18) and (19)), involve these functions for all relative times 
from — oo to +oo. Thus these equations are extremely difficult to handle in practice. A 
drastic approximation is to put the relative time equal to zero, i.e. to write an equal-time 
formalism. This route has been followed in other approaches 0. This is equivalent to the 
so-called quasi-particle approximation, in which the energy dependence of the spectral 
function is approximated by a 5 function, and thus off-shell properties of the distribution 
functions are neglected. As mentioned in the introduction this is a drastic approximation 
for many reasons. An improvement is therefore expected by expanding all quantities 
around r = 0. The off-shell information is then contained in the time derivatives with 
respect to r at r = 0, either inside of the convergence radius by Taylor expansion or 
by analytic continuation. We therefore will develop an infinite set of equations, or a 



hierarchy of equations, for these equal-time derivatives originating from the kinetic and 
the constraint equations. 

The equal-time r-derivatives of order n g^"'\T, r) are related to the energy moments 
of order n of the Fourier transformed function g{T, u) since we can write 

^(")(T) := d:g{T,T)\r=o 

= ^—^Jdcucu-giT,u;). (21) 

However, the energy moments may not exist in general. E.g. if g{T.,uj) behaves like 
-^ for large u , then the second and higher energy moments of the spectral component 
diverge. Such a behaviour e.g. occurs in the gradient approximation of the spectral 
function neglecting the dependence of the imaginary part of the retarded self-energy on 
uo. The reason is that integration and differentiation cannot, in general, be interchanged 
in eq. (21). In contrast to the energy moments the time derivatives g^"'\T) exist as seen 
from the definition of g in terms of operator products. 

The problem of the divergence of the energy moments could be overcome by in- 
troducing a cut-off in the energy-integral, but we do not choose this procedure here. 
However, we can retain the physical intuition that, e.g. for the spectral Green func- 
tion, the zeroth derivative g^^\T) is related to the normalization of the function, the 
first derivative g^^\T) to the mean energy of the quasi-particle, and the second deriva- 
tive g^'^\T) to the width of the spectral function. Indeed, neglecting the possible di- 
vergence problems, we have for the spectral function A defined e.g. by Henning [|l| 
A{T,uj) = -^Jdre'^^g'^iT-^r). Then, / rftucjM(T, w) = i"+i^^W and the n*'^ r- 
derivative is exactly the n^^ energy moment of the four dimensional Wigner transform 
of the spectral function up to a factor i""*"^. 

We start by setting up the hierarchy for the spectral function g"^. From the upper 
row of eqs. (20) for g^ we obtain the hierarchy of kinetic equations for the equal-time 
r-derivatives of the spectral component. The equal-time commutation relation, through 
eq.(9), results into 

/^°^ = -^ (22) 

which is the sum rule for the spectral function. Using this relation, the first three 
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equations of the kinetic hierarchy are 

^^5^(°)-2^sini^(e-/(°))=0 

.^/(^) - 2^ sin \0{e ■ /«) - cos \o{^ ■ /(«)) = 

^^/(^) - 2^ sin l^(e ■ /(^)) - 2 cos \<){^ ■ /«) = 

2zsin^^(a^(0)./(i))_,^ (23) 

2 

where e = ^ + ^"5 (T) and where we used the identities 

P-V^/ = -2sini^(^-/) 
m 2 2m 

-(^Ar-p^)/ = -2 cosier—-/) 

The integrals on the rhs of eqs. (20) for g^ vanish at r = 0, eqs. (18). /,Froni the second 
row of equation (20) we obtain a constraint hierarchy from which we write down the first 
two members. 

-2cosi^(e-(7^(°)) + 22(7^« = 

-2 cos -^(e ■ (?^(i)) + 2^(7^(2) = 2 cos io(a^(°) ■ (?^(°)) (24) 

In general the n*'* equation of the constraint hierarchy gives the n*'* r-derivative 
of the spectral function, (yi^^^"-), in terms of the lower order r-derivatives of the spectral 
function, g'^^'-''' and of the spectral self-energy , a^^^^ with j < n — 1. On the other hand 
the (n + 1)*^ equation of the kinetic hierarchy is a kinetic equation for the n*^ r-derivative 
of the spectral function, g^^^\ in which only the lower order r-derivatives of the spectral 
function and spectral self-energy, g"^^^' and a"^^^' for j < n — 1, appear. The first kinetic 
equation (23) is satisfied by (yi^C^) = — i, eq. (22). 

The natural question arises: are the n^^ equation of the constraint hierarchy and 
the (n + ly^ equation of the kinetic hierarchy compatible? In fact, the (n + I)*'* kinetic 
equation is not a new equation, but it can be expressed as a combination of lower order 
constraint and kinetic equations. Since we have not made any approximation yet, this 
result is true beyond the mean field approximation, in which only the singular part of the 
self-energy, i. e. Hartree-Fock self-energy, is retained. We will show this explicitely in the 
following up to order n = 2. This means for the spectral part that the algebraic solutions 
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of the constraint hierarchy satisfy the kinetic hierarchy, and that the time evolution is 
described only by the T-dependence of the self-energies. 

With eq. (22) and eq. (24) we have from the constraint hierarchy 

^^(1) = -e (25) 

<7^(2) = ^cosi^(e-e)-(7^(°) (26) 

To prove the compatibility of the kinetic and constraint hierarchies we have to 
verify that the solutions (25,26) satisfy eqs. (23). This is done in Appendix B using 
the properties of the differential operator ^ given in Appendix A. Therefore, only the 
constraint hierarchy is necessary, the kinetic hierarchy being redundant, as the solutions 
of the former satisfy automatically the latter one. The r-derivatives of the spectral 
function can be determined from the constraint hierarchy ij the r-derivatives of the 
spectral self-energy are known. In general, these depend on the r-derivatives of the 
kinetic part of the Green function, g^ ^ as will be shown later . Thus the spectral and 
kinetic Green functions are coupled. 

4 Hierarchies of equations for the kinetic Green func- 
tions 

Following the same way as for the spectral function we obtain from eqs. (20) two infinite 
hierarchies of equations for the equal-time r-derivatives of the kinetic part g^ . Now also 
the equal-time derivatives of the collision integral enter. The first three equations of the 
kinetic hierarchy are 

,^^^(0) _ 2z sin i^(e ■ g""^'^) = e^^^i^^^^ - e-^^Sf^^^ 
t^g^W - 1i sin -^(e ■ (?^(i)) - cos -^(1^ ■ (?^(°)) = 



^^g^(^) - 2i sin -0{e ■ g^^^^) - 2 cos -<}{^ ■ g^^'^) - - sin -0{^ ■ g""^'^) = 
dT 2 ^ ' 2 ^dT ' 2 2 ^dT^ 



where S^ are given by eq. (19). 
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The first two equations of the constraint hierarchy are 

-2 cos -0(e ■ (7^(°)) + 2^<?^(i) = e^^^f ^^^ + 6-^^55^°^ 

_ , sin 1^(^ . ^^(0)) _ 2 cos ^Oie ■ g^^'^) + 2tg^^'^ = e'^<^S'^^'^ + e'^^^^^^^ (28) 

Let us first consider the mean field approximation, when colhsions are neglected 
{S^{t) = 0). The n*^ equation of the constraint hierarchy gives the n*'^ r-derivative, 
g^m^ in terms of the lower order r-derivatives, g^^-^' for j < n — 1. The {n+iy^ equation 
of the kinetic hierarchy is an evolution equation for the n*'* r-derivative, g^^"'\ In this 
case we show in Appendix B that at least up to order n = 2 only the first equation of 
the kinetic hierarchy and all the equations of the constraint hierarchy are independent; 
the n^^ equation of the kinetic hierarchy, n > 1, can be expressed in terms of the lower 
order kinetic and constraint equations. Therefore the two hierarchies are compatible in 
the mean field approximation (at least up to order n = 2). We note that in general the 
singular part of the self-energy on the contour may depend on p, in contrast to the usual 
assumption 0]. 

The inclusion of the collision contribution is more difficult because we have to 
approximate the integrals S± and S± . The main contribution to the integrals 

5f (°) = ± f_^ dt[a^{T + I Tt)g''{T + | ±t) - a^'(T + | ±t)/(T + | ^t)] 

comes from small values of the variable of integration t since large values involve Green 
functions at large relative time, which are expected to be small. Therefore we can 
approximate the integrals using as a lower limit a judiciously choosen parameter —A 
instead of — oo (or Tq — T). For this purpose we expand the functions under the integral 
of the type f(T + 1/2, ±t) in t around f(T,t = 0), i.e. around the equal-time value at 
current gobal time T, and integrate over dt. If we expand the integrand to first order in t 
we obtain the second order approximation in the parameter A for the collision integrals 

+ -A2(a^(^)^^W - a^(o)^^(i) - a^(i)/W + a^(o)/(i))] , (29) 

where the time argument T has been omitted everywhere. For the derivatives of 5*:^^ and 
5*5, by expanding the integrand in zero order in t, we obtain an expression linear in A 
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^Ml)gK{0) ^ ^A(0)gK(l) _ ^K{l)gA{0) _ ^K{0)gA{l) 



2\ dT ^ dT dT ^ dT J 



(30) 



Another way to approximate the colhsion integrals by a convergence factor e~" is pre- 
sented in the appendix C There it is seen that the results are essentially the same as 
above, with A replaced by a term proportional to 1/a. 

Using the above expressions (29) and (30) for the collision intregrals the first two 
equations of the kinetic and constraint hierarchies for the kinetic Green function, g^ , 
include memory effects due to the presence of time derivatives. The eqs. (29) and (30) 
are second, respectively first, order polynomial approximations in A for the collision 
terms. In general, in the /c*^ order approximation of S*.]- in A appear r-derivatives of 
the Green functions and of the self-energy up to order k — 1. As we will see in the next 
section, the latter can be obtained in terms of only the same order r-derivatives of the 
Green functions. However, in the same approximation in A, the collision terms in the 
n*^-order equations of the two hierarchies contain r-derivatives of the self-energy up to 
order n + k — 2. Thus the appropriate truncation of the hierarchies is not as obvious as 
in the mean field approximation; only the hierarchy for the spectral function truncates 
by itself, completely similar to the mean field case. 

If we approximate the collision terms S^ in order k, k > 1, with respect to A 
then the first equation of the kinetic hierarchy contains the r-derivatives of the Green 
functions up to order k — 1. Therefore in general we choose to approximate the collision 
terms which appear in the n*^-order equations of the two hierarchies, S± , in order 
k — {n—l). Within this approximation the first k — 1 equations of the constraint hierarchy 
contain also the r-derivatives of the Green functions up to order k — 1. Therefore, one 
has k equations for k r-derivatives. The r-derivatives of the self-energies appear up to 
order k — 1 , i.e. N r-derivatives and they can be obtained self-consistently. 

In the following we will consider the k = 2 approximation. The kinetic and con- 
straint equations are 

zA^^(°) - 2z sin 1^(6. ^^(°)) = 

2A [l - ^-pi cosi^(a^W^^(°) - a^(°)^^W) 



12 



+zA2 sin -0 (a^(i)(?^(o) - a^(o)(^^(i) - a^(^)(?^(o) + a^(o)(?-^(i)) (31) 

2zA (l - ^A^ sini^(a-^(^^(°) - a^(%^(°)) 

+A2 cos -0 p(i) (^^(°) - a-^(o)(?^(i) - a^^'^g^^'^ + ^^^(^-^(i)) (32) 

The familiar Vlasov and Boltzmann equations are particular cases of the above 
approximation. When in the collision integrals S^ only the terms that are zeroth order 
in A are retained (the collision integrals are neglected), the equations for g^^^> has the 
same form as the corresponding equation for ^^(°) . Thus the first equation of the kinetic 
hierarchy is 

^^^^(°)-2.sini^(e-^^(°)) = 

which in the gradient approximation is the Vlasov equation with a self consistent mean- 
field, which is included in e. The validity of this approximation is not directly related 
to the quasi-particle approximation. Here we only use the assumption that the collision 
term can be approximated in zeroth order in A. 

The lowest order nontrivial approximation is A; = 1, which results in taking only 
terms with g^^^^ and g^'^'^^ in eqs. (31,32). We then have the kinetic equation for g^^^") 

,^^^(0) _2^ sin i^(e- (7^(0)) = 2Acos io(a^(°)(7^W - a^^V^^)) . (33) 

Let us note that the term (^a^^^'g^^^' _0-^{o)(^^(o)^ jj^ ^^le above equation can be put into 
the familiar form of the collision term of a Boltzmann-like equation, i. e. 

^^A(0)gK(0) _ ^K(0)gA(0)^ ^ (^>(0) _ ^<(0))(^>(0) + ^<(0)) _ (^>(0) ^ ^<(0))(^>(0) _ ^<(0)) 

= 2(a>(°)(?<(°) - a<(°)(?>(°)) (34) 

Eq. (33) goes beyond to the gradient approximation, and the higher order space deriva- 
tives contain information about quantum effects. In the gradient approximation we 
obtain a Boltzmann-like equation for the zeroth r-derivative of the Green function. 

Consider now the k=2 approximation for the collision integrals, eqs. (31, 32), we 
need to calculate the first two r-derivatives of the kinetic Green function. In the gradient 
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approximation (cos ^<) -^ 1, sin |<) — > |{, }) eq. (32) is an expression for g^^^^ depend- 
ing on g^^^^ and the first order r-derivative of the self-energy. This expression can be 
substituted into the kinetic equation (31) yielding a transport equation for g^^'^^ Since 
g^^^^ appears in eq. (31) under the Poisson bracket, the Poisson bracket in eq. (32) has 
to be neglected consistently with the gradient approximation, because a substitution of 
it would result in a term in eq. (31) which involves second order derivatives. Moreover, 
g^^^^ appears in eq. (31) in the second order term in A. Therefore we can use in eq. (31) 
the solution of the algebraic eq. (32) for g^^^^ in zeroth order in A g^^^^ ~ —ieg^^'^\ 
Thus, in a consistent approximation, substituting this expression in eq. (31), we obtain 



where cr^^^-* can be calculated using the same approximation in for g^^^^ in terms of g^'^^^ 
as above. 

To go beyond the gradient approximation one has to keep the higher order gradient 
terms in the otherwise algebraic constraint equation for g^^^\ 



5 Self-energy in the Born approximation 

To solve the equations for the r-derivatives of the Green functions we need the corre- 
sponding r-derivatives of the self-energies. The self-consistent Hartee-Fock self-energy is 
singular in time on the contour; it appears in the equations as a potential as (momen- 
tum dependent in general) and it depends only on the zeroth r-derivative of the Green 
functions. The Hartree self-energy can be written |12| 



S,^(l,l') = -25(ti-ti05(xi-xiOy'rf3a;2l^(xi-X2)G<(x2,ti,X2,ti) 

The Wigner transformation does not depend on momentum due to the presence of the 
^-function in space coordinates and it also has no r-dependence. One obtains 

aH(R,T) = -zJd'^V{Oj^Jd'pg^^'\K-^,p,T) 

' Y.—,dR\ld'P9^^'\^^V,T)]fd'irV{i) (36) 



(2vr)= 
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The Green function appearing above is 

g<'\R,p,T) = 



W 



^W(R,p,T)-l 



(37) 



using definitions (7) and the sum rule (22) for the spectral function. The gradient 
expansion in the second line of eq. (36) gives non-locality contributions to the Hartree 
self-energy due to the dependence of the Green function g^^^^ on R which are not present 
in a translationally invariant system. However, the lowest non-locality contribution 
(proportional to the gradient of the Green function) has a factor / d^C,C,V{C,) which is zero 
for inter-particle interaction potentials with V{^) = V{—C,)- Therefore, the specification 
of the Hartree self-energy in the local approximation 



(7h{'R,T) 



i2ny 



d'pg<'\R,p,T) 



d'^ViO 



(38) 



is justified at least in the gradient approximation. The second order term in the gradient 
expansion in eq. (36) for spherically symmetric potentials V{^) is proportional to 



A 



R 



d^g<'^\R,p,T)\Jd'^eV{0. 



(39) 



For a (5- interaction potential the local approximation is exact; as the nucleon - nucleon 
potential is of short range, one can expect that in nuclear processes the local approxi- 
mation is reasonable even for large non-localities. 



In the same manner we obtain for the Fock self-energy |T2 

i 



d^^e-'P^ViO I d^p'g<°^ (R, p', T)e-'P'^ 
^^^^Jd'p'g<'\R,p',T)vip'-p), 



(40) 



where v is the Fourier transform of the interaction potential. The above expresion has 
no gradient corrections. Thus, the Fock contribution is the convolution of the Green 
function g"^^^^ with the Fourier transform of the interaction potential. We note that the 
Hartree- Fock self-energy can be obtained from only zero order r-derivatives of the Green 
function g^^^^ using eq. (37). 

The lowest order contribution to the non-singular part of the self-energy on the 



Schwinger-Keldysh contour is given by the direct and exchange Born diagrams |jT2|. We 
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will discuss here only the contribution for the direct term, the exchange term is given in 
an analogous way. 

a>f{K, p,T;T) = jlyJ d^e-'^^ J d'p, J d' P2V {p{)V {p^) j dW"""'^ 
g>'<iR, p^, T; r) j d'p2e'^-^^-^^^^P-^^g>^<{n + ^i^, p^, T; r) 
/dV-^P^^«-(^^+''^»^<'>(R+ ^^,P3,T; -r) (41) 

Expanding eq. (41) around r = 0, it is seen that the n^^ r-derivative of the Born self- 
energy involves the r-derivatives of the Green functions up to order n. This property 
remains true if one replaces the Born approximation for the self-energy by the T-matrix 
approximation if the off-shell values of the T-matrix can be approximated by the on-shell 



ones as it was discussed in |jT3|. Further we make a gradient expansion in eq. (41) of the 
last two Green functions around R. In the local approximation, i.e. in zero order, we 
have for the zeroth r-derivative of the direct Born self-energy 

a|,<(°HR,p,T) = ^/rf3p2/rf=^P3^>'<(°nR,P2,T)^<'>W(R,P3,T) 

^>'<(°)(R, p - P2 + P3, T)v\p2 - P3) (42) 

In a choosen order in the gradient expansion for the constraint and kinetic equa- 
tions we need the self-energy in the same order. It is straightforward to prove that the 
first order terms in the gradient expansion of eq. (41) for r = cancel. Thus, as for the 
Hartree-Fock part of the self-energy, the local approximation for the Born self-energy 
is consistent with the gradient approximation for the transport equations. When con- 
sidering second order terms in the gradient expansion for the transport equations, also 
second order terms in the expansion of eq. (41) have to be retained. The contribution of 
second order gradient terms to the zero order r-derivative of the direct Born self-energy 
is of the form 

3>-<"\n, p - P2 + P3.r)aft5«2(p2 - p^) (43) 

In the local approximation, from eq. (42) with eqs. (7) and eq. (22) one has 
.a^r (R, P.T) = \-^J d'p, J d'p, { [^^(°)(R, P2, T)^^(°)(R, pa, T) + i; 
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+ <7^(°)(R,p-p2 + P3,T)[^^W(R,P2,T)-(7^W(R,P3,T)]}t;2(p2-p3) 



1 1 



J d'p2 J d^ps [l + ^^(°)(R, P2, T)g^^'\R, P3, T)] t;2(p2 - Ps) (44) 



4 (27r)6 

In the second step we assumed that g^^^\p) = g^^^\—p) which is vahd e.g. in the 
fireball region of a heavy ion collision. This step is not necessary but we use it for the 
following estimate. Using the non-interacting Green functions [|T| 

^>W = _,(l_^p)) , ^<(°) = m(p) (45) 

where n{p) is the occupation number of the state with momentum p and < n(p) < 1 
we obtain an estimate of icr^^'^'' in eq. (44) 

n(P3)ll-n(P2)]]t''(P2-ps) (46) 

We note that the spectral part of the self-energy in the Born approximation is an integral 
of a positive function and thus positive. 

Eqs. (25,26) in the gradient approximation (cos(i<) -^ 1) imply 

g^{r) = -le-'^--'-^-' (47) 

where we approximate ln5f'^(r) in second order around r = 0. Here the spectral function 
decays by a gaussian in the relative time when the spectral part of the self-energy is 
calculated using eq. (46). Then a characteristic correlation time is 



1 
n 



ZCTq 



e 

This estimate can be used to fix the parameter A (or a in Appendix C) in the calculation 
of the collision integrals, eqs. (29,30), choosing A = ri (or a = y —^ — ) where the 
spectral self-energy, ia^ , is determined as above using the initial occupation numbers. 

6 Numerical results 

We used eq.(33) (k=l approximation of the collision integrals) to test the influence of 
the terms not usually retained in the gradient approximation on the time evolution of 
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the phase space nucleoli density distribution, n = (1 — ig^^^')/2. To keep the numerical 
problem as simple as possible we restrict ourselves to a one-dimensional configuration 
space (the phase space is two dimensional); physically, this could be the case of two 
colliding nuclear matter slabs. For simplicity, a nucleon-nucleon interaction of the shape 
of an attractive rectangular well potential with typical depth (40 MeV ) and range (2 
fm) was used. 

We discretize the partial differential equation, eq. (33), and solve the resulting 
finite difference equation. We evolve the phase space density distribution in time steps 
of 0.1 fm/c ; the density distribution is given at mesh points between -8 fm and +8 
fm in steps of 0.4 fm in space and between -8 fm"^ and +8 fm~^ in steps of 0.4 fm~^ 
in momentum. The initial density was taken as a superposition of two gaussians (the 
widths were 2 fm in position and 1 fm"^ in momentum) centered at -2 fm and +2 
fm in position and at +2 fm^^ and -2 fm^^ in momentum; the mean momenta of the 
gaussians correspond to an energy of the colliding slabs of 80 MeV/ A. 

The contour plot of the results of various approximations are collected in Fig. 
1. The initial phase-space density is presented in panel (a). To test the program we 
calculated first the phase-space density for free motion, i.e. when there is no interaction 
and only the first order gradient terms are retained from the drift term in the left hand 
side of eq. (33) and the collision terms in the right hand side of eq. (33) are neglected. 
The results after 10 fm/c are presented in panel (b) and they are consistent with the 
classical Hamiltonian fiow {i.e. with the test particle method). 

The next step is to include the mean field (the Hartree-Fock self-energy). The 
results obtained at 10 fm/c in the gradient approximation (Vlasov equation) are pre- 
sented in panel (c) of Fig.l. It is seen that there is considerable distortion of the phase 
space distribution. Including the next order gradient correction in the Vlasov equation 
the results are given in panel (d). We find that the quantum corrections to the drift 
term result in a much smoother behaviour of the phase-space distribution function, in 
fact, more similar to the free evolution. The results of including the collision terms but 
with the first order gradient expansion is given in panel (e). The hole in the center of 
phase space is filled and the distribution is more extended. This calculation corresponds 
to the usual BUU approach. In panel (f) both the next order gradient and the collisions 
terms are included. The density distribution becomes still smoother. 
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It is seen that higher order gradient corrections result in a smoother phase space 
distribution. This might seem contradictory since with higher order gradient corrections 
one moves closer to the evolution of a quantum state, the Wigner transform of which is 
expected to show more oscillations. However, it has to be remembered that the initial 
distribution does not correspond to a pure quantum state. 

In Fig. 2 we present the corresponding density distribution in space and mo- 
mentum; they are obtained by integrating the phase-space distribution over momenta, 
respectively over space. One can see that the spreading in space of the density distri- 
bution increases by including the third order gradient term in the drift part (panel (d) 
in Fig. 1) relative to the Vlasov approximation (c) and it increases further when the 
collisions are included in the calculation (e,f). The momentum distributions are rather 
stable with respect to the gradient terms (d), but the collision terms result in a consid- 
erable spreading of the distribution in momenta (e,f), as expected. With a closer look 
one sees that in the Vlasov approximation the momentum distribution spreads in a self 
consistent mean-field (c), but this is no longer the case when the gradient (quantum) 
corrections are taken into account (d). 

These calculations include nonlocal terms simultaneously in the drift and the col- 
lision part of the Bolzmann-like equation (33). The inclusion of nonlocality terms in the 
scattering integral given by Morawetz et al. [^] leads to consistency with the thermody- 
namic virial corrections, but we think that nonlocal corrections should be also included 
in the drift term for the case of nonequilibrium systems. 

Up to now we have discussed higher order gradient terms, i.e. quantum corrections 
to the space evolution. It is of interest to look also at the off-shell properties of the Green 
functions, i.e. at the properties of the spectral function. A first look is obtained by 
discussing g'^'^'^\ i.e. the width of the spectral function. According to eq. (26) it depends 
on (T"^(°) which we have already calculated. According to eq. (48) it is also proportional 
to the inverse of the parameter A, connected with memory effects which are discussed 
below. The width of the spectral function is shown for the slab-slab collision in Fig. 3 as a 
function of the one-dimensional distance variable for different times, in particular for the 
initial time (panel (a) of Fig. 1) and for 5 and 10 fm/c for the BUU approximation (panel 
(e) of Fig. 1). It is seen that the width increases in the time evolution of the initially 
nonequilibrated momentum configuration. It is largest during the phase of maximum 
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overlap of the slabs. In the later phase of the evolution it becomes more spread out and 
lower. This is an expected behaviour as the final configuration should evolve towards hot 
nuclear matter (the initial collective kinetic energy transforms in thermal energy) and 
the final width has a contribution from thermal effects. It is to be noted that the width 
goes to zero automatically outside of the system because there the spectral function has 
to go back on shell. 

Another quantum effect on the time evolution of density distribution is due to 
the inclusion of higher order terms in A in the collision integrals. This is a memory 
effect as it involves time derivatives of the self energies. To illustrate numerically this 
quantum effect we used eq. (35) in the simplest case, namely two interpenetrating nuclear 
matter systems in one dimension, i.e. relative to the previous case the space extension 
is disregarded. In such a translationally invariant system all the quantities depend only 
on momentum and time and eq. (35) can be written for the the time dependent density 
in momentum space n 

^^^ = 2A(p;T) (l - ^^^) [(l-n(p;T))(-^a<(°)(p;T))-n(p;T)(.a>(°)(p;T))] 

(49) 
where we used the relations g^^^^ = 2i{n — 1/2) and (y(^(°) = — 1 and the self-energies in 
the Born approximation are 



^a>(°)(p;T) = .^|rf3p2/rf'p3t^'(P2-P3)[l-^(p2)]n(p3)[l-n(p-p2 + P3)] 
-2a<(°)(p;T) = -^y'rf'p2/rf'p3^^'(P2-P3)n(p2)[l-n(p3)]n(p-p2 + p3)] 



and, using eq. (48) 



^^P' ^^ " V(.a>(0)(p;T)) + (-.a<(0)(p;T)) ^^°^ 

If we neglect the second order term in A we obtain an equation which has the structure 
of a Boltzmann-like equation 

^ = IgainiV; T) - kossiV, T) (51) 

with the collision integrals of the Uehling-Uhlenbeck form 

Igainiv; ^) = J ^^^2 j d^PsS^P, P2 " Ps; T) [1 - n{p)] [1 - n{p3)]n{p2)n{p - P2 + P3) 
Iiossip; T) = J d^p2 J d^P3s{p, P2 - Ps; T)n{p) [1 - n(p2)]n(p3) [1 - n(p - P2 + P3)] 
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but with an effective in- medium cross-section s(p, P2 — Pa; T) that depends not only on 
the transferred momentum, but also on the actual distribution. Thus the above equation 
includes nonequilibrium effects through the in-medium cross-section 

s = -(W^j d'^P'2 J d%v^{p2 - P3)[n(p3)(l - n(p2)) + n(p - p'2 + P3){n{p^) - ^(pg))] 

which is a complicated functional of n{p;T). Thus nonequilibrium effects affect the 
collision integrals, in the sense that the cross section is dependent on the momentum 
distribution. In general it is therefore not justified to use a cross-section calculated 
for equilibrated configurations in momentum space in a BUU description of a nuclear 
collision. 

When we retain the second order term in A in eq. (49), using eq.(50), we obtain 



IgainiP; T) - IlossiP] T) (52) 



2 lI^(-^'°'(p^^»-4(Ij 

In equilibrium the rhs of eq. (49) vanishes and one obtains 

2 

'A(^T) 



.a<(°)(p;T) = n(p;T) (53) 



2 



<^°np;T) = (i-Mp;^))x(^ (54) 

Let us remark that if we assume that the condition (53) remains true also in nonequi- 
librium (this can be justified if the system is not too far from equilibrium) and use it in 
the Ihs of eq.(52) we obtain exactly eq. (51). More generallly we can write 

2 

'A(^T)2 

where C is a correction term of the form 



la 



<(o). 



p;T) = n(p;T)— — - + C(p;T) (55) 



C(p;T) = .^—^ J d^p^ J d^p3v\p2 - Ps) {n{p2)[l - n{ps)]n{p - P2 + P3)]- 

n(p; T)[n(p3)(l - ^(pa)) + n(p - P2 + P3)(^(P2) - ^(pg))]} (56) 

and eq.(52) can be rewritten as 

^ = Igain{P; T) - IlossiP] T) ^ fff^"^' ' ^ ' 

which differs from eq.(51) only in the last term which can be considered as a memory 
term. 
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We investigated numerically eq.(51) and (57) in a one- dimensional case starting 
with an initial distribution of two gaussians. In the case n{p) = n{—p) the parameter 
A does not depend on p, but it depends on time. The results are presented in Fig. 4. 
One can see that using A calculated at the initial time instead of a dynamical A in 
eq. (51), i.e. without non-equilibrium effects on the in medium cross-section, the high 
momentum distribution tails are overestimated. Further inclusion of the memory term, 
eq. (57), changes mainly the distribuiton at small momenta. 

7 Conclusions 

A transport theory for equal-time quantities related to the spectral and the kinetic part 
of the one-particle Green function has been derived. The present description is not 
in principle restricted to the mean field approximation and higher order terms in the 
gradient expansion can be systematically retained. 

We note that in the mean field approximation the above derived equations reduce 
to the full quantum equations of motion for the Wigner function in a self-consistently 
determined mean-field. We can then study the effect of higher order terms in the gradient 
expansion, i. e. the quantum corrections to the classical equations of motion. This could 
be of interest for systems which are localized in space, as is the case for heavy ion 
coUissions, especially in the initial stage. The condition of validity for the gradient 
approximation, the product of the characteristic lengths at which the Green function 



varies in position and in momentum is much larger than 1 [T^, is not fulfilled for a 
nucleus (one is of the order of the nuclear dimension, ~ 1 fm, and the other of order of 
Fermi momentum, ~ 1 frrT^). 

We also consider the inclusion of terms which go beyond the mean-field approxi- 
mation in the self-energy and make a detailed discussion of the self-energy in the self- 
consistent Born approximation. The compatibility of the two hierarchies (kinetic and 
constraint) for the r-derivatives of the spectral Green function at zero relative time 
(equivalent to the energy moments of the spectral function provided they exist) has 
been proved up to the second order. This extends the similar result obtained in the 
mean-field approximation in [^]. 

For the kinetic part of the Green function similar hierarchies of equations have been 
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derived and the possibility of their systematic truncation is discussed. The lowest order 
truncation is a Boltzmann-like equation (33) for the equal-time kinetic part of the Green 
function; all the terms in the gradient expansion are retained and this allows to study, as 
in the mean field case, the effect of the terms neglected in the gradient approximation. 
The collision term in eq. (33) has the Uehling-Uhlenbeck form (see eqs. (42) and (34)), 
in which the Pauli exclusion principle is automaticaly taken into account. 

Let us mention that it is possible to prove that any truncation of the gradient 
expansion in eq. (33) conserves the total number of particles if the first order one, 
i.e. the usual gradient approximation, does. Integrating eq. (33) over position and 
momentum, and using integration by parts to transfer all the derivatives to the Green 
functions (all the selfenergies are supposed to vanish together with all their derivatives 
at infinity, and A is considered constant) the only surviving term in the rhs of eq. (33) 
will be 2AJd^RJd^p{a^^^^g^^^^ — a^'^^^g^^^^), i.e. the term present in the gradient 
approximation. In the gradient approximation, the change in the Boltzmann equation 
in the next order truncation of the hierarchies has been derived, eq. (35). 

The methods develloped in this work were tested in a simple model of one-dimensional 
slab-on-slab collisions. Here we went only to the lowest order of the hierarchy, i.e. we 
remained in the quasi-particle approximation. However, we investigated the effect of the 
next order gradient terms, i.e. of true quantum corrections, and found them of consid- 
erable importance. Thus these corrections could also influence the dynamical evolution 
of real heavy ion collisions. 

The next step will be to consider higher order truncations of the hierarchies. This 
will result in a dependence of the time evolution of phase-space density on the higher 
order moments of the Green function g^ . As these moments are a measure of the off- 
shellness, this could be an alternative way to include the off-shell motion in heavy ion 
collision calculations. For other proposed approaches in the literature see |^, p. 

The approach of the present paper can be used to treat a system of bosons, or 



a coupled fermion-boson system, as quantum hadrodynamics [|T6| which is the modern 
way for describing the nuclear many-body problem in the framework of effective field 
theories. 
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Appendix A 

Properties of the differential operators sin |<) 

and cos ^(} 

e^^ U ■ 9) = f -*^ g 

is a noncommutative, associative product defined on tlie space of complex valued func- 
tions /(R, p). 

It is necessary to verify the associativity only for functions of the form fi = 
gi(aiR+ftp) y^[^\^ arbitrary ctj, Pi as every function can be expressed as a linear com- 
bination of the above functions ( Fourier decomposition) with complex coefficients and 
the product -k is obviously distributive. Using 

it is a straightforward exercise to prove that 

/i * (/2 * /s) = (/i * /2) ^ /s • 



We have the identities 



1 1 

cos-^(/-^) = -{f-kg + g-kf) 

1 1 

sin -<){f ■ g) = — (/ -kg-g-kf) 

and the following relation is obvious 

sin^^(/-/) = [Al) 

Using the definitions of the differential operators and the associativity we obtain 

sin -<>(/ ■ cos -^{g ■ h)) = cos -^^g ■ sin -<}{/ ■ h)) + cos -<>(/i • sin -<>(/ ■ g)) 
The particular case of the above relation g = f is 

cos ^<>if ■ sin ^^(/ ■ h)) = sin ^^(/ ■ cos ^^(/ ■ h)) {A2) 

We have also the identity 

11 11 11 

cos -<}{g ■ cos -<>(/ ■ h)) + sin -<}{g ■ sin -<>(/ ■ h)) - sin -<>(/ ■ sin -<}{g ■ h)) 
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= ^[9-^{f-^h + h-kf) + {f-kh + h-kf)-kg-g-k{f-kh-h-kf) 

+ {f'kh-h-kf)-kg + f-k{g-kh-h-kg)-{g-kh-h-kg)-kf] 

= ^[f-^{9-^h + h-kg) + {g-kh + h-kg)-kf] 



cos^<)(f . cos^O(g . h)) {A3) 



Appendix B 

Compatibility of the constraint and kinetic hierarchies for the Green 
functions 

Using eq. (22), eq. (Al) for / = e, ^ = ^ and cos i^(/ ■ ^^(o)) = fg^^^^ = -if 
it is a straightforward algebraic exercise to prove that the second kinetic equation (23) 
is satisfied by the solution (25) of the first constraint equation. 

To show that the solution of the second constraint equation (26) satisfies the third 
kinetic equation (23) we take the action of the operator i-^ — 2isin|^(e-) on the eq. 
(26) and we obtain 

^^/^^^ - 2^ sin l^(e . /(^)) = -2 cos i^(§ ■ e) - ^^ 
+2 sin i^(e ■ cos ^0{e ■ e)) + 2t sin ^<}{e ■ a^(°)) 

= -2 cos ^OQ . e) - ^^ - 2^ sin ^0{a^i^^ ■ e) , ^^^^ 

where in the last step we used the identity (A2) in the Appendix A 

sin-^(e-cos-<>(e-e)) = cos -<>(e ■ sin -<>(e ■ e)) = 

The third equation (23) can be rewritten using the relation (25) 



^^/^^^-2.sinl^(e./(^)) = -2cosi^(§.e)-.^ 

-2isin-^(a^(o)-e) 



(52) 



Comparing eqs. (Bl) and (B2) we conclude that the solution (26) of the constraint 
hierarchy satisfies the third equation of the kinetic hierarchy. 
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In the same way we will show the compatibility of the constraint and kinetic hi- 
erarchies for the kinetic part of the Green function in the mean field approximation. 
The second kinetic equation (27) can be obtained from the first kinetic equation and the 
first constraint equation (28). Acting with the operator i-^ — 2isin|<)(e-) on the first 
constraint equation we obtain {S^ = 0) 

^^^^«-2^sini^(e-^^«)+^(^^-2^sini^(e))cos^^(e.^^(°)) = 

Using the first kinetic equation we obtain 

^^,^« - 2^ sin i^(e ■ ,^«) - cos i^(§ ■ ,^(°)) - 
2cosi^(e-sini^(e-(?^(°))) + 2sini^(e-cosi^(e-(?^(°))) = 

The last two terms cancel and we obtain exactly the second kinetic equation (27). 

In the same way one can obtain the third kinetic equation. Acting with the operator 
2isin |<^(e-) on the second constraint equation (28) and using the first two kinetic 

equations we obtain 

^-^g^^'^ - 2^ sin l^(e ■ ,^(^)) - cos lo(§ ■ ,^(^)) - 

i cos ^0{e ■ cos i^(^ ■ ^^W)) + 2 sin ^0{e ■ cos ^0{e ■ (?^«)) + 
^sini^(e.sini^(^.^^W))-^sini^(^.sini^(e.^^(°))) = 

The 5*^ and the 7*^^ cancel and we obtain 

. A,^(2) _ 2^ sin l^(e . ,^(^)) - '- sin i^(|^ • ,^(°)) - cos 1^(§ • g^^^^) = 
-^ jcos i0(6 . COS ^-0{^ ■ ^^W)) + Sin ^-0{e ■ sin ^-0{^ ■ ^^W))- 
sin i^(^ ■ sin i^(e ■ ^^(°))) | = -^ cos ^^(^ • cos i^(e ■ ^^(°))) 

where we used the identity (A3) in Appendix A. Using the first constraint equation (28) 
we obtain 

,_|_^^(2) _ 2^ sin -^(e ■ g^^'^) - - sin -Oi^ ■ g""^'^) - cos -<>{% ■ g^^'^) = 
dT 2 ^ ' 2 2 ^dT^ 2 ^dT ' 
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which coincides with the third kinetic equation (27). 



Appendix C 

Approximation of integrals appearing in the hierarchies for g^ 

In the hierarchies of equations for g^ the terms due to collisions, S± are of the 
form / = J_^dxf{x). We will approximate such an integral introducing an arbitrary 
real and positive parameter a 

rO 



dx 



e-^^f\x) 



^ax 



If we assume that e~"^/(a;) is a function localized around zero (this is reasonable for the 
functions discussed here depending on the relative time) and a is choosen such that the 
radius of convergence of the Taylor series for / around zero is greater than - {a > -) 
then the main contribution to the integral comes from the region in which the function 
is well approximated by its Taylor series in the origin and we have 

1 <9" , __ „, .. , /-o 

ra=0 

Using the notation f^^'^ = -^f{x) |^=o we obtain 



Z^n n\ ox"- J-oo 



1 oo n 1 r(k) 

EE ^ 



We consider truncations of the sum over n, Y^n=Qi ^^'i "^^ "^ill ^se A^ = 1 or A^ = 2 to 
approximate the integral by the derivatives of / up to order N . 

We note that if the radius of convergence of the Taylor series is infinity then the 
sum of the above series gives the exact result which does not depend on a. However, the 
result of the truncation of the series depends on the value of the parameter a; a low order 
approximation is good if - is choosen as a characteristic decay time of /, physically a 
correlation time of the system. For example, if we calculate in the N = 2 approximation 
the integral with / = e~^ ^ using a = P the approximation will be 4 which differs from 
the exact value ^4 by 10 percent. 
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We can use the above results to evaluate the integrals Sj^. and S^ . It is 
straightforward but tedious to write the collision terms in order N = 2. The results are 

a 2a oT 

and 



1 
± — 

3a 



(^^A(l)gK(0) ^ ^A(0)gK{l) _ ^K{l)gA{0) _ ^KiO)gA{l)-^ 



2V (9T ^ (9T dT ^ dT J 

The similarity of these expansions with the A-expansion used in Chapter 4 is 
obvious. It corresponds to replacing - by A with slightly different numerical factors. 
The reason is also clear: in the A-expansion we cut the integral over the past history 
with a sharp cut-off A, in the present method by a smooth exponential cut-off with a 
parameter a. 
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Figure 1: Contour plots of the phase space distribution of slab-on-slab coUisions at 
80 MeV per nucleon displayed for the direction of the finite extension of the slabs (x- 
direction). Panel (a) shows the initial distribution. The other panels show the distri- 
butions at time 10 fm/c in the following approximations: (b) free motion; (c) Vlasov 
approximation, i.e. first order gradient expansion; (d) Vlasov plus next (third) order 
gradient terms; (e) Vlasov plus collision term (corresponds to the usual BUU approx- 
imation); (f) Vlasov plus next order gradient plus collision terms. The contour lines 
are give in panels (a - d) for 0.2(0.2) and in panels (e,f) for 0.1(0.1) with the notation: 
outermost contour line (increment). 
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Figure 2: Space (top) and momentum (bottom) distributions for the same slab-on-slab 
collisions as in Fig. 1, i.e. the projections on the coordinate and momentum axes. The 
curves correspond to the following approximations (with the appropriate panel of Fig.l): 
free motion (solid, panel b); Vlasov (longer dashed, panel c), Vlasov plus next order 
gradient terms (short dashed, panel d), Vlasov plus collision terms (dotted, panel e), 
Vlasov plus next order gradient plus collision terms (dash-dotted, panel f). 
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Figure 3: Width (i.e. second order energy moment) of the spectral function at different 
times as a function of the coordinate of the one-dimensional model. 
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Figure 4: Time evolution of an initial distribution of two gaussians in momentum 
(top) and at 10 fm/c (bottom). The curves corresponds to the following approxima- 
tions: Boltzmann with time dependent (nonequilibrium) cross-section s(p, P2 — Ps; T) 
(solid line), Boltzmann with time independent cross-section calculated from the initial 
configuration s(p,p2 — P3',T = 0) (long dashed line), Boltzmann with time dependent 
cross-section plus memory term (short dashed line). 
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